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We calculate exactly the von Neumann and topological entropies of the toric code as a function of system 
size and temperature. We do so for systems with infinite energy scale separation between magnetic and electric 
excitations, so that the magnetic closed loop structure is fully preserved while the electric loop structure is 
tampered with by thermally excited electric charges. We find that the entanglement entropy is a singular function 
of temperature and system size, and that the limit of zero temperature and the limit of infinite system size do not 
commute. The two orders of limit differ by a term that does not depend on the size of the boundary between the 
partitions of the system, but instead depends on the topology of the bipartition. From the entanglement entropy 
we obtain the topological entropy, which is shown to drop to half its zero-temperature value for any infinitesimal 
temperature in the thermodynamic limit, and remains constant as the temperature is further increased. Such 
discontinuous behavior is replaced by a smooth decreasing function in finite-size systems. If the separation of 
energy scales in the system is large but finite, we argue that our results hold at small enough temperature and 
finite system size, and a second drop in the topological entropy should occur as the temperature is raised so 
as to disrupt the magnetic loop structure by allowing the appearance of free magnetic charges. We discuss the 
scaling of these entropies as a function of system size, and how the quantum topological entropy is shaved off 
in this two-step process as a function of temperature and system size. We interpret our results as an indication 
that the underlying magnetic and electric closed loop structures contribute equally to the topological entropy 
(and therefore to the topological order) in the system. Since each loop structure per se is a classical object, 
we interpret the quantum topological order in our system as arising from the ability of the two structures to be 
superimposed and appear simultaneously. 



L INTRODUCTION 

Some strongly correlated quantum systems have rather rich 
spectral properties, such as ground state degeneracies that are 
not related to symmetries, but instead to topology. 1,2 Such sys- 
tems are said to be topologically ordered, 3 and they can have 
excitations with fractionalized quantum numbers, 4 as in the 
case of the fractional quantum Hall states. There have been 
proposals to utilize topologically ordered states for fault tol- 
erant quantum computation, exploiting the resilience of these 
systems to decoherence by local perturbations or disturbances 
by the environment. 

Levin and Wen 5 , and Kitaev and Preskill 6 recently pro- 
posed that a characteristic signature of topological order can 
be found in a subleading correction of the Von Neumann (en- 
tanglement) entropy in systems prepared in (one of) its ground 
state(s). This topological correction to the entanglement en- 
tropy was indeed confirmed by exact calculations in discrete 
models exhibiting topological order, as well as in contin- 
uum systems such as fermionic Laughlin states. 7 The notion 
of topological entropy provides a "non-local order parame- 
ter" for topologically ordered systems. Hereafter, we refer to 
topological order as characteristically identified by such non- 
vanishing topological entropy. 

Although quantum topological order was introduced as a 
pure zero-temperature concept, is was recently shown 8 that a 
closely related behavior can be observed also in mixed state 
density matrices that describe classical systems in the pres- 
ence of hard constraints. These findings show that topolog- 
ical order can survive thermal mixing under certain condi- 
tions, e.g., in hard constrained systems. Moreover, any pos- 
sible experimental observation of quantum topological order 



must take into account the fact that the T = limit is only 
an idealization and temperature, albeit small, is a perturba- 
tion that cannot be neglected. This is particularly relevant, 
for example, if one is interested in a practical application of 
topological order towards quantum computing, which will al- 
ways be done at finite temperature. It is therefore interesting 
to study the behavior of topologically ordered systems as the 
temperature is gradually raised from zero, in search of a uni- 
fied picture of topological order encompassing both the quan- 
tum zero-temperature limit and the classical hard-constrained 
limit. 

In this paper, we investigate the fate of quantum topological 
order in the two-dimensional toric code on the square lattice 
in thermal equilibrium with a bath at finite temperature. In 
particular, we do so by studying the entanglement and topo- 
logical entropies of the system, which we compute exactly. 

We start from the zero-temperature limit of the model, 
which has been thoroughly studied in Ref. 9. In this limit, the 
ground state (GS) of the system can be mapped onto two loop 
structures 10 each of which, we argue, is responsible for half of 
the topological contribution to the von Neumann entropy (i.e., 
half of the topological entropy of the system). As temperature 
is raised from zero, thermal equilibration disrupts (breaks) the 
loop structure and it is expected to destroy topological order. 
With an exact calculation in the limit where one of the two 
loop structures is fully preserved while the other is allowed 
to thermalizefAA we show that the topological entropy gradu- 
ally decreases as a function of temperature, for fixed and finite 
system size, from its zero-temperature value down to precisely 
half of that value. In particular, the temperature dependence 
of the topological entropy can be shown to appear always 
through the product K A (T) N, where K A {T) is a monotonic 
function of temperature with K A (Q) = and K A (oo) = oo, 
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and N is an extensive quantity that scales linearly with the 
number of degrees of freedom in the system. Therefore, the 
thermodynamic limit TV — > oo and the T — > limit c/o nof 
commute, and if the former is taken first, the topological en- 
tropy becomes a singular function at T = 0, and it equals one 
half of its zero-temperature value for any T ^ 0. In other 
words, in the thermodynamic limit any infinitesimal temper- 
ature is able to fully disrupt any loop structure for which we 
allow thermalization, and the contribution from this structure 
to the topological entropy is completely lost (irrespective of 
the presence of a finite energy gap). On the other hand, finite 
size systems can retain a statistical contribution to the topo- 
logical entropy (in the sense that its value varies continuously 
with temperature) originating from a thermalized underlying 
loop structure. 

From our results, we then infer the behavior of the finite- 
temperature topological entropy in the generic case, as illus- 
trated in Fig. Q] For finite size systems, we expect to ob- 
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FIG. 1 : Qualitative behavior of the topological entropy as a function 
of temperature T and number of degrees of freedom 27V, for the 
generic case where the two coupling constants in the model are well 
separated, i.e., X A <C X B . The exact shape of the first crossover is 
shown in Fig. [5] 

serve two continuous decays of the topological entropy, due 
to the gradual disruption of each of its two loop contributions. 
Each drop occurs when the number of corresponding defects 
~ N e~ 2X A,B/ T reaches a value of order one, where 2N is 
the number of degrees of freedom in the system, and X A and 
X B are the two coupling constants in the model, associated 
with one loop structure each. The separation between the two 
decays (the quantity A in the figure) is therefore proportional 
to the difference between the two coupling constants X A and 
X B in the Hamiltonian. Once again, if the thermodynamic 
limit is taken first, both decays collapse into a singular behav- 
ior where the topological entropy vanishes everywhere except 
for T = 0, where its value depends on the order between the 
thermodynamic and zero-temperature limit. Notice also that, 

although the crossover temperature Tcmm ~ b /\xlVN 
goes to zero in the limit of N — > oo, it does so only in a 



logarithmic fashion. 

The case of classical topological order is recovered in the 
present study when thermal fluctuations are allowed to com- 
pletely break one of the loop structures while the other is 
strictly preserved. Indeed, this can be accomplished by impos- 
ing appropriate (local) hard constraints on the classical analog 
of the toric code. 8 Our results illustrate how both the concept 
of quantum topological order and of classical topological or- 
der are equally fragile, as they truly exist only in the zero- 
temperature / hard-constraint limit. Their effects however can 
extend well into the finite-temperature / soft-constraint realm 
- as our calculations show - so long as the size of the system 
is finite. 

Our results suggest a simple pictorial interpretation of 
quantum topological order, at least for systems where there 
is an easy identification of loop structures as in the case here 
studied. The picture is that (i) the two loop structures con- 
tribute equally and independently to the topological order at 
zero temperature; (ii) each loop structure per se is a classical 
(non-local) object carrying a contribution of In D to the topo- 
logical entropy (D = 2 being the so-called quantum dimen- 
sion of the system); and (iii) the quantum nature of the zero- 
temperature system resides in the fact that two independent 
loop structures can be superimposed (therefore leading to an 
overall topological entropy equal to 2\nD = hiD 2 ). In this 
sense, our results lead to an interpretation of quantum topo- 
logical order, at least for systems with simple loop structures, 
as the quantum mechanical version of a classical topological 
order (given by each individual loop structure). 

We also investigate the von Neumann (entanglement) en- 
tropy S VN as a function of temperature and system size. For 
instance, we show that, given any bipartition (A, B) of the 
whole system S = A U B, the quantity 
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where Sy N (T) is the entropy of partition A after tracing out 
partition B, m B is the number of disconnected regions in B, 
and L = y/~N is the linear size of the system. From this result, 
we learn that the topological contribution to the entanglement 
entropy can be filtered out directly from a single bipartition, 
provided m B > 1, as opposed to the constructions in Refs.ljlg 
that require a linear combination over multiple bipartitions. 

We also show that, as soon as the temperature is different 
from zero, the von Neumann entropy is no longer symmet- 
ric upon exchange of subsystem A and subsystem B, and it 
acquires a term that is extensive in the number of degrees of 
freedom that have not been traced out (see Eq. d39ll). Symme- 
try and dependence only on the boundary degrees of freedom, 
at least in the thermodynamic limit, can be recovered if one 
considers instead the mutual information 
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as we explicitly show in this paper. Notice that I A g{0) = 
^vn(O) = 'S'vn(O)- Once again, the mutual informa- 
tion exhibits a singular behavior at zero temperature since 
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the thermodynamic limit and the zero-temperature limit 
do not commute. We find that the explicit topologi- 
cal contribution to I AB (T) in the thermodynamic limit is 
— h (m A + m B — 1) In 2, where m A (m B ) is the number of 
disconnected components of partition A (B). 

Although we consider here a very specific model, we be- 
lieve that our results are of relevance to a broader context, at 
least at a qualitative level. For example, it would be interest- 
ing to investigate the specific behavior of systems where the 
underlying structures responsible for the presence of topologi- 
cal order are no longer identical to each other. This is the case 
of the three-dimensional extension of the toric code, where 
a closed loop structure becomes dual to a closed membrane 
structure.— 

The paper is organized as follows. In Section|Il]we present 
the model and discuss its characteristic features and proper- 
ties. In Section [Til] we compute the von Neumann entropy of 
the system as a function of temperature and system size, in 
the limit of one of the coupling constants going to infinity. 
We then obtain the exact expression for the topological en- 
tropy in Sec. [IVj and we illustrate its behavior with numerical 
results. Finally, we discuss the implications of our results for 
the system with finite coupling constants and we infer the full 
temperature and system size dependence of the topological 
entropy in Sec. [V] Conclusions are drawn in Sec. [VU 



II. THE FINITE-TEMPERATURE TORIC CODE 

The zero-temperature limit of the model considered here 
was studied by Kitaev in Ref. |9l It can be represented by 
2 A spin- 1/2 degrees of freedom on the bonds of an L x L 
square lattice, N — L 2 , with periodic boundary conditions 
(toric geometry). The system is endowed with a Hamiltonian 
that can be written in terms of star and plaquette operators as 



H = -A, 
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E A * 
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plaquettes p 



where A^ and X B are two positive coupling constants, B p = 
Yliszp °t> an d A s = Iljgs °j> wilh i labeling all four edges of 
plaquette p and j labeling all four bonds meeting at vertex s of 
the square lattice. Notice that the Hamiltonian, all the B p op- 
erators and all the A s commute with each other, and one can 
diagonalize them simultaneously. Given that there are A — 1 
independent plaquette operators and A — 1 independent star 
operators (ITj,=i B p = I — Ilfli A B ), the eigenvectors with 
fixed B p and A s quantum numbers form a 2 2 -dimensional 
space. Furthermore, one can show that the GS 4-fold degen- 
eracy has a topological nature that can be split only by the 
action of non-local (system spanning) operators. The ground 
state wavefunctions of this model are known exactly, 9 and can 
be written in the cr z basis as 



l*o> = 



1 



3|0), 
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where G is the Abelian group generated by all star operators 



the dimension of G, and |0) = \a\ ... er z ) is any given state 
that satisfies the condition B p | 0) = |0),Vp. There are four in- 
equivalent choices for |0), corresponding to the four different 
topological sectors of the model. The choice of sector is im- 
material to the results presented hereafter, since they all have 
the same entanglement, 13 and we will set |0) = | + + . . . +) 
for convenience throughout the rest of the paper. 

Notice that the system is symmetric upon exchange of <7 X 
with a z components and of stars with plaquettes on the lat- 
tice. In the cr x -basis, each site must have 0, 2 or 4 spins with 
a negative er x component on the adjacent bonds. If we were 
to remove all the bonds with a negative er x component, we 
would obtain a configuration of closed loops on the square 
lattice, where loops are allowed to cross but do not overlap. 
Once a convention is established on how to interpret sites en- 
tirely surrounded by spins with positive <r x component (e.g., 
as two different loop parts touching at the corner, say the up- 
right and down-left loops), then one can establish a one-to-one 
correspondence between all basis states and all loop configu- 
rations on the square lattice where loops cannot overlap and 
can at most touch at a corner in an up-right, down-left fashion 
(see Fig. |2). The same is true for the a z -basis, but the loops 
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FIG. 2: Illustration of the spin-loop correspondence discussed in the 
text. All vertices with two positive and two negative <r x components 
on the adjacent links can be obtained via appropriate rotations of the 
ones shown in the figure. 

now live on the dual lattice given by the centers of the plaque- 
ttes of the original lattice. This description of the GS of the 
toric code in terms of loop degrees of freedom gives a qualita- 
tive picture of the origin of the non-local behavior of the sys- 
tem from which the presence of topological order stems. In 
particular, given that the two simultaneous loop descriptions 
can be mapped one onto the other upon exchanging cr x with 
<7 Z as well as the square lattice with its dual, it is tempting to 
speculate that they contribute equally to the topological order 
present in the system, and each loop structure is responsible 
for precisely half of the topological entropy. Our exact calcu- 
lations show that this naive picture is indeed correct: if either 
of the loop structures is exactly preserved while the other is 
destroyed, e.g., via coupling to a thermal bath, the topological 
entropy of the system lowers to half of its original value. 
Based on the er x and a z loop description, all possible per- 
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turbations to the system can be qualitatively divided into three 
different classes: (i) those that couple to a crMike term, and 
are able - if sufficiently strong - to disrupt the underlying a x 
loop structure, but not the er z one; (ii) those that couple to 
a tr x -like term, with precisely the opposite effect; and (iii) 
those that couple to a cr y -like term, and are able - again, if 
the coupling constant is large enough - to disrupt both loop 
structures, thus leading to a vanishing topological entropy. A 
generic coupling to a thermal bath is likely to encompass all 
of the above terms and in the thermodynamic limit the vanish- 
ing of the topological entropy is unavoidable. For finite size 
systems and at low enough temperature, however, the relative 
scale of the two coupling constants X B and plays a cru- 
cial role in determining how effective each of the above terms 
is with respect to the others. In this paper we consider the 
case when the two energy scales are well separated, namely 
X B 3> A^, and we discuss qualitatively the behavior of the 
system as the separation becomes weaker and vanishes. A 
large separation between the two energy scales is indeed ex- 
pected if we notice that the toric code is a lattice realization 
of a Z 2 gauge theory, where the two coupling constants A^ 
and X B relate directly to the chemical potential of the electric 
and magnetic monopoles. 14 On the ground of a large separa- 
tion between the two energy scales in the Hamiltonian, three 
distinct temperature regimes can be outlined: 

(a) T -C X A J In y/N, when all thermal excitations have 
a small Boltzmann weight and for finite size systems 
at finite time scales the topological entropy effectively 
retains its zero temperature value because of the scarcity 
of defects that can disrupt the loop structure; 

(b) X A /]ny/N < T < X B /\n\TR, when thermal excita- 
tions of the <x z type can disrupt the <r x loops structure, 
while the er x -like excitations are rather unlikely to occur 
and they can be effectively neglected; and 

(c) X B / In \/~N < T, when the appearance of all the three 
types of thermal excitations leads to the complete dis- 
ruption of the topological contribution to the entangle- 
ment entropy. 

(Notice that the opposite case, where X A X B , leads to 
equivalent results based on the symmetry of the model.) The 
temperature range considered in this paper corresponds to 
regimes (a) and (b), where the a 1 loop structure is effectively 
preserved for sufficiently small system sizes and time scales. 
We will then discuss how our results can be used to infer the 
behavior of the topological entropy across the whole temper- 
ature range, illustrated in Fig.[T] 

Basically, one can define a temperature dependent defect 
separation length scales £ AB ~ e A - 4 -- B ^ T ', so that as long as 
the sample size is below these corresponding scales, the sys- 
tem is free from the associated type of defects. Similarly, one 
can define temperature-dependent time scales for defects to 
appear. The toric code is fragile in the sense that 0(1) defects 
destroy its topological order, so that for practical considera- 
tions not only the temperature must be small compared to a 
gap, but the system size and the time scales must not be too 
large as well. 



We now consider the simplification where the finite system 
length \/N <C £ b , so we can neglect defects in the cr z loop 
structure. Forbidding any defects in the <r z loop structure is 
equivalent to neglecting all thermal processes that violate the 
constraint Yiiep a \ = +1' ^P- Therefore, the Hilbert space in 
the regime of interest and within the chosen topological sector 
(recall that |0) = | + + ...+» is given by {g\0) | g e G}. 
The equilibrium properties of the system are then captured by 
the finite-temperature density matrix 

P{T) = \ e-& 

E 9eG (0|ffe-^. 9 |0) ' ( } 

where we used the group property to write a generic element 
g" 6 Gas g" = gg', 3\ g' S G given g g G. Recall that all 
group elements, as well as their composition, are defined as 
products of star operators modulo the identity Yl^=i A s = 

For the model under consideration, it is convenient to 
rewrite the Hamiltonian as 

H = -X B P-X A S 

plaquettes p 

S = ]T A s . (6) 

stars s 

Notice that Pg\0) = Ng\0), Vfif € G, and therefore 

(0\ge-P H gg'\0) = e^* N (0\ge^ s gg'\0) 

= e (3X * N (0\e fsx A s g / \0) (7) 

where we used the fact that any g commutes with S by con- 
struction. 

Now, recall the definition of a group element g E G, which 
can be symbolically represented by the notation 9 = X\ seg A s 

(modulo the identity n«=i A s = 1, i.e., g = Yl s< z g A s — 
Yisdg Ag). Given the expansion 

e px A S = -Q [ cosh/3AA + sinh/JA^ A s ] , (8) 

stars s 

which follows from the definition S = 5Z s =i A s and from the 
fact that A 2 S = 1, one obtains 

(Ole^VlO) = 

= (0| J[ [cosh/3A A + sinn^ A,] ]\ A s , |0) 

s s'eg' 

= (cosh f3X A ) N (tanh/3A j4 )" (9 ' ) 

+ {cosh/3X A ) N (t&nhpX A ) N - n(9,) , (9) 

where TV is the total number of stars in the system, and n(g') 
is the number of flipped stars in g' . Notice that the ambiguity 
in the definition of g = Y[ seg A s = Yisfg A s ~ namely the 
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fact that if g' is given by the product of a set of A s , it is also 
given by the product of all other A a but for those in the set 
- does not affect the equation above. In fact, this ambiguity 
amounts to the mapping n(g') <-> N — n(g'). Similarly, 



Z = ]T(0| ff e- 

geG 



g\o) 

n[ 
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Substituting Eq. (0 and Eq. ( fTOb into Eq. (O after relabel- 
ing K A = — ln[tanh(/3A yl )] gives 



P(T)= E 



1 



<g') + e -K A (N-n(g')) 



\G\ 



[1 



-K A m 



9 



99 



s.s'eG 

(ID 

In the limit of T — > (/3 — > oo), if^ — * + , all g' are equally 
weighed, and one recovers the density matrix of the zero- 
temperature Kitaev model. In the limit T — > oo ((3 — * 0), 
if^ — > oo, all <?' are exponentially suppressed except for 
g' = 1, and one recovers the mixed-state density matrix of 
the topologically ordered classical system discussed in Ref.H. 
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where we used the generic tensor decomposition |0) = \Q A ) (8 

I°b). 9 = 9a® 9b- and the fact that (°b\9b9b9b\°b) = 1 if 
g' B = 1 B and zero otherwise. The latter follows immediately 

from the fact that the group G is Abelian and that A 2 S — 1, V s, 

and therefore 3g = 1 g for any choice of £?. We also denoted 

by G A = {g € G | g e = 1 B } the subgroup of G given by 

all operations g that act trivially on B (similarly for G B in the 

following). For convenience of notation we defined 
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III. THE VON NEUMANN ENTROPY 

Let us consider a generic bipartition of the system S into 
subsystem A and subsystem B (S = A U B). Let us also 
define (S g ) to be the number of star operators A s that 
act solely on spins in A (B), and T^ AB as the number of star 
operators acting simultaneously on both subsystems. Clearly 
these quantities satisfy the relationship Y, A + E g + ^ab ~ 
N. Whenever a partition is made up of multiple connected 
components, e.g., A = A 1 U ... U A m with A i n Aj =0 
and A i connected, Vi, j, let us denote with Y, A the number 
of star operators acting solely on A i (Ea = J2i ^a )• [Since 
in the following we will consider only the case where either A 
or B have multiple connected components, but not both at the 
same time, will be used unambiguously to denote the 

number of star operators acting simultaneously on A { and on 
B or on B i and on A, according to the specific case (S^ B = 

Y.,>-AB '-I 

The von Neumann (entanglement) entropy 5 VN of a bipar- 
tition (A, B) is given by 



5y N = -Tr [p A \np A ] = Sy N , 



(12) 



where p A = Ti B (p) is the reduced density matrix obtained 
from the full density matrix p by tracing out the degrees of 
freedom of subsystem B, and the last equality holds whenever 
the full density matrix p is a pure-state density matrix. 

In order to compute the von Neumann entropy ( flZb from 
the finite-temperature density matrix (fTTT i. we first obtain the 
reduced density matrix of the system using the same approach 



Notice that a star operator A s can either act solely on 
spins in partition A (represented in the following by the no- 
tation s E A), solely on spins in partition B (s G B), or 
simultaneously on spins belonging to A and B (which we 
will refer to as boundary star operators, and represent by 
s € AB). As discussed in Ref. 8, a complete set of gener- 
ators for the subgroup G . can be constructed by taking: (i) 
all star operators that act solely on A, i.e., {A s \ s € A}, 
together with (ii) the collective operators defined as the prod- 
uct of all stars acting solely on a connected component of 
B times the product of all boundary stars of that specific 
component, for all the m B connected components of B, i.e., 
{IlseB. x Ils'e.4B. ^s'j V connected components i}. No- 
tice that not all the collective operators are new operators with 
respect to those generated by the star operators in A. In fact, 

IL (iL'efi, A s> Us"eAB t A s») = UseA A s' and one can 
show that there are precisely m B — 1 new, independent op- 
erators. Consequently, the cardinality of the subgroup G A 
is given by d A = \G A \ = 2 s ^ +m s~ 1 . Similarly for G B , 

d B = \ G b\ = 2 S B+''n A -\ 

To proceed with the calculation of the von Neumann en- 
tropy of the finite temperature system, it is useful to use the 
above set of generators in order to represent the group G A 

/- T . ■ -ii c n r\ ~\ connected components i 

in terms of Ising spin variables {0 S , &i} seA , 
where 9 S = —1 (1) corresponds to the star operator A s ap- 
pearing (not appearing) in the decomposition of g G G A , and 
similarly Q i = —1 (1) corresponds to the collective opera- 
tor n se B ! A s x Us'eAB l A s> appearing (not appearing) in 
the same decomposition. Notice that the correspondence is 
2-to-l, since a configuration {9 S , s and its spin-flipped 
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counterpart {8 S , 0j}j ;S , where 9 S = —6 S and 0^ = -0^ 
map onto the exact same g e G A (which follows from the 
fact that one of the collective operators can be generated out 
of the others appropriately combined with the star operators 
in A). 

In this representation, 
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and the restricted summation S/gTo e (0j is subject to the 
constraint (O^^ A . . . g' n _^\0 A ) ^ 0, which can be explicitly 
stated in terms of the spins 8 S and 0^ as 
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t.S / — 1 



i,5 /— 1 



1=1 



E^ + JE 5 ^ 6 ^), d5b) 



TV 1 
T + 2 



where we used the fact that + J2i ypB + ^,48 J = N 
and we introduced the notation Y, v = E e + ^ab ■ 

Let us then use Eq. (13[ to compute the n-th power of 
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9i,a\ a)( a\9i,a9i,a92,a\ a)(°a\92,a9'2,a 

■■■9n,A\ A)(°A\9n,A9n,A- ( 16 ) 

Each expectation value above imposes g l A g[ A = 1^, 

/ = 1, ... , n — 1, and therefore g l g[g l+l £ G B . Upon re- 
labeling n — 1 summation variables so that = <7;<?/5/ +1 
for / = 1, . . . , n — 1, the corresponding sums can then be 
combined with the respective inner product and they can be 
written as J^. eG 1 = d B . Therefore, the equation above 

can be simplified to 



Pa(T) 



E II E •fe-cso 



9i,a\°a)( a\9i,a3i,A ■ ■ ■ 9n,A- ( 17 ) 



Taking the trace of p A (T), using the fact that all the g's 
commute, and J2 g ggl = |G|, one obtains 



w\) n e >?t(sd<o^u---<a> 

, v n—1 1 n constr. 

E *r({«M«>}), (is) 



where the factor of 1/2™ comes from the 2-to-l nature of the 
representation of G A in terms of Ising spin configurations, 



Above and in the following, the short-hand notation i, s in 
sums and products stands for connected components i and s S 
A. 

Let us then substitute Eqs. ( fT~5T > into Eq. ( fT4t . 



co S h(^E s ^ + ^E,^ 4 0f) 



cosh 



2 cosh ( ^A 



E< 

j=±i 



2 cosh 



e n« 

J=±i se^i 



4jEx, e (!) 



1 , (19) 



and take the sum over all possible {9s , 0^ } i s configura- 
tions (without any constraint), 



2 cosh ( ^-N 



e e n 

j=±i {6/< !) ,e (!) }. seA 

e fn e 

/ 



e 2 



e 2 



n E e^ J ^ e '° 



(20) 



ey=±i 



Using the expression above, one can rewrite Eq.[T8las 



7 



B 

|G| 



2 2, 



cosh 



HE HE 



n e 



n 5 di e " } - i)^n w - 1) + n <n ^ + *. m + 



i,s i= 1 



Ls l—l 



1=1 



and, upon expanding the product n[Li j =±i ^(^z )) = 2{j }p Il/Li )> we °btain 



B 

|G| 



2 2, 



cosh 



E 



( 

n e 

V n 



e 2 



■E, J,t 



J 



n e e_ 

4 {e! !, }r=i 



to 



n E « 

V ne< !, =-i 



n e e 

1 {e<°}r =1 



/ J 



(21) 



Notice that the indices s and i to the variables Os and 9^, re- where Zn /a> (K A , { J l }) is the partition function of a chain of 
spectively, are mute since the sums over the possible values of n Ising spins with periodic / antiperiodic boundary conditions, 
9 { s l) and are performed first. Therefore one can simplify in presence of a nearest-neighbor interaction with position- 
the notation above by replacing # -> 0(0 and 0« - 9« dependent reduced coupling constant X A J, /2. Similarly, 
for all s and z. Moreover, one can recognize that an Ising spin 
chain {9^} with Y[ ^ — ±1 is dual to an Ising spin chain 
{t ; } with periodic / antiperiodic boundary conditions (using 
the 2-to-l mapping 6"> = r ; t ;+1 ). Thus, 



(p/°) / 



_^ p. /a. 

^ ^ g~2 4 "El = l ''l^' ' — _ ^ ' g — 2^ El = l ^1 T ! r l + l 



E 

{e«} r=1 
ne< ! )=±i 



(23) 



{ e(o }?=1 
ne< ! >=±i 



= -^/ a) (iT A ,{JJ), (22) and Eq. (EB can be rewritten ; 



2 2 ™ | cosh ( 



E >< 



2 S ^+™ B \\G\J 2 2n L 



(24) 



The partition functions can be evaluated, for either bound- 
ary conditions, using a transfer matrix approach. For periodic 



boundary conditions, we obtain 



Z^(K,{J l })=Tr 



.i=i 



T, 



e 2 J i e 2 J i 

_i£ J K J 

g 2 I g 2 i 
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Since all matrices T ; are diagonalized by the same unitary 
matrix 



For the antiperiodic case, 



1/11 



and using the cyclic properties of the trace, we get 



(25) 4 a )(^,{Jj)=Tr 



Z^(K,{J t }) = Tr 



.1=1 



Tr 



Tr 



11 / K 



n 



e 2 1 + e 2 



e 2 J < — e 2 1 



-q ( 2cosh(§) 



i=i 



2 J, sinh(f ) 



2 cosh 



K 



n I n 



2 sinh 



A" 



"iII T « 



i=i 



1 
10/' 



and we get 



Z^iKAJt}) = Tr 



K 

2 cosh ( — 



1=1 

n 

in 



i=i 

n / n 



\l=l 



2 sinh — 
2 



(27) 



(26) We can substitute into Eq. d24l) and we obtain 



Tr[#] 



cosh 



cosh 



d B 2 E ^+ r 



2 
K 



yln 



1=1 / 

n 



cosh (^w) 
mi)]) ( 



sinh 



E >< 

n| 
n 



cosh 



cosh 



A, 



A, 



n / n 



n 

ip 



sinh 



sinh 



\i=i 



2 

A, 



Notice that the factors Y\i Ji can be dropped since the two terms between curly brackets in the equation above get simply 
exchanged by the two different values of J l = ±1. The summation over { J l }f =1 gives therefore just a multiplicative factor 
2 n . Recalling the definition d A = 2 s ^ +m e~ 1 , we finally get to 



Trht/V 



d A d B 



2cosh(^fiv) ycosh(^iv) 



|G| 



cosh 



cosh 



A, 



sinh 



A, 



A, 



sinh 



A 



n 



cosh 



cosh 



sinh 



A, 



A. 



sinh 



A 



(28) 



We can now use the identity 



simplified notation: 



lim3 n Tr[#] 

n— >1 



(29) 



cosh 



K A 



sinh 



K, 



(30a) 



to compute the von Neumann entropy Sy N (T) = 
— Tr [p A \np A ]. It is convenient to introduce the following 




(30b) 
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which allows us to rewrite the terms in the last two lines of and 
Eq. ([28]> as 



i 



(31) 
(32) 



d n F ± 



(n) 



One can verify that 



^ A e ± ^ {N - 1) (x\nx±y\ny) 

(34) 



-I- -* /v 

e 2 



(33) 



Using this notation 



d A d B 



■2>nsh(^N) ycosh^V) 



|G| 



}■ 



and 



2 cosh 



In 



1 d^dg 



cosh 



|G| 



2 cosh ( —^-N 



+11^^^-^ (xlnx + y\ny) + J2 e * 



d A d 



= - ln -j^p + In cosh ( -^-AT 



E^e^^" 1 ) (xlnx + ylny) 



2 cosh (^iV*) 

-^e-^ 1(A '- 1) (xlnx-t / lnt / ) + ^ -V^-E 



hiXi + y, In&) 



ln^_£ +lncosh ( -^JV) - E^ (xlnx) 



cosh(^f(iV- 1) 



cosh 



sinh(^-(7V- 1) 



cosh 



cosh(^(iV- E p )) sinli(^(iV-Sp 



cosh 



cosh 



(35) 



For a finite size system, the limit of K A — > + yields In the limit of iT^ — > oo instead, x, y ~ e KA ^ 2 /2, x i , y 



x, x , — > 1, y, yj — > and therefore 



S4(T^0)^-lnV B 



|G| ' 



(36) 



consistently with the known zero-temperature result 
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10 



, K A s r i / 2 j2 t anc [ therefore 



5<? N (T^oo) 



, n d A d B 

\G\ 



K, 



■JV-ln2-S 



A 



K, 



-In 2 



Et, -ln2 



ln^-ln2 

|G|' 



E^ ln2- 



m e In 2 



(37) 



where = 2 s ^ +m s _1 . This result is consistent with Ref.|H 
where the classical system that obtains in the T — ► oo limit 
had been previously discussed. 

Let us then consider the thermodynamic limit L — ► oo. The 
presence of terms that depend on the product between K A and 
extensive quantities such as N, E ,,, and E p requires care- 
ful consideration. As we will see, they will indeed lead to 
a singular behavior of the von Neumann entropy, as well as 
the topological entropy discussed in Sec. [IV] For any finite 
K A £ (0, oo), the limit L — ► oo (N = L 2 — ► oo) yields 



^vnCO 



|G| 



^A-ln2 



-E_^ (xlnx + j/lny) e 2 



(38) 



We can further simplify this expression in the limit 
Ep ,S, 3> 1, i.e., in the limit of large partitions; 15 this is 
the case of interest, for example, in the definition of the topo- 
logical entropy discussed in the Section [IV] If we expand 

x i) y i ~ e <AT " p ^ 2 /2, we obtain 



&vn (O 



■ d A d B 

|G| 



K, 



N - In 2 



-E^e 2 (a; In a; + y In y) 



E 



K, 



In ^r 4 ^ 6 + (m £ 



In 2 



l)ln 2 



-E, 



|G| 

e 2 (imi + j/m?/) — 



= (^AB -™a) ln2 

— E_4 e 2 (a; ln a; + y In y) - 

which is consistent with the limit Sy N (T 



K, 



(39) 



— ln(dg/|G|), but no longer consistent with the finite size re- 
sult 0) = - ln^dg/IGQ. In fact, 



lira SUT) = -ln%^ + (m B -l)ln2, 



(40) 



while 



lim SUT) = -ln%^. 

L^oo,T^0 |G| 



E^g - m A — m B + 1 ) ln 2 



(41) 



Notice that already in Eq. ( 1381 the limit T — ► fails to give 
the known result, Eq. ( |36l l. However, it is only in the assump- 
tion that all partitions scale with the system size that we arrive 
at Eq. ( |39l . Notice also that the difference between the two 
orders of limits arises only if subsystem B has more than one 
connected component. Moreover, such difference is of order 
one, while the common term — \n{d A d B )/\G\ scales linearly 
with the size of the boundary between A and B (oc E^ g ). 
Therefore, from a thermodynamic point of view the order of 
limits is immaterial to the von Neumann entanglement en- 
tropy. On the other hand, we will see in the following section 
how this difference plays a crucial role in the topological en- 
tropy of the system, and gives rise to a singular behavior at 
zero temperature in the thermodynamic limit. 



A. The mutual information 

In the previous section, and in particular from Eq. ( f35T > we 
clearly see that the von Neumann entropy is symmetric upon 
exchange of A and B, and it satisfies the so-called area law 
only at zero temperature (and for topologically ordered sys- 
tems, symmetry is lost unless the T — ► limit is taken before 
the L — ► oo limit). At any finite temperature, Sy N (T) acquires 
an extensive contribution scaling with the number of degrees 
of freedom in partition A. In this sense, the von Neumann en- 
tropy ceases to be a good measure of the entropy contained in 
the boundary between the two subsystems, and therefore is no 
longer a good measure of entanglement at finite temperature. 

Alternatively, one can consider another quantity called the 
mutual information I AB , which we redefine for convenience 
multiplied by a factor 1/2, 



Iab(T) 



1 



^vn(0 



"5vn(0 



qAUB 



CO 



(42) 



oo 



so that I AB (0,N) = 54(0, AT) = 5y N (0, N). Substituting 
the result from Eq. d35l >. and recalling that d AuB = \G\, d* — 
1 and E_4 UB = N, we obtain the behavior of I AB (T) as a 
function of temperature and system size for a topologically 
ordered system, 
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The mutual information has the immediate advantage over the 
von Neumann entropy that it is symmetric upon exchange of 
subsystem A with B. Moreover, while the general expression 
above seems to have an explicit dependence on the bulk de- 
grees of freedom (Ei and S e ), we show hereafter that in 
the thermodynamic limit as well as in the zero-temperature 
limit the mutual information depends only on the boundary 
degrees of freedom (S^ e ) and, as such, is a good candidate 
for a measure of entanglement at finite as well as zero temper- 
ature. 

In the T — > limit we recover the known result 



Wo) 



. i n d -* d B 



\G\ 



= £ 



J AB 



1 In 2. 



(44) 



In the T — > oo limit instead we obtain 



1 In d -* d B 

2 |G| 



-'AB ~ m A 



1 In 2, (45) 



which again shows only a dependence on boundary degrees of 
freedom, with the addition of topological terms of order one, 
and is symmetric under the exchange of A and B. 

We can then compare these finite-size results with the be- 
havior of 7_ AZ3 (T) in the thermodynamic limit L — > oo for 
large partitions, 

d A dr. 



Iab(T) 



In 



L A U B 



\G\ 



-'AB 



4(- 



J AB 



'A 



(x In x + y In y) e 
+ m K - 1 ) In 2 



K 



4(- 



2 In 2 + (x In X + y In y) e - — 



',4 



1 B 



1 In 2. 



Once again this is consistent with the finite size limit T 
but not with the T — > limit. In fact, 



(46) 

> oo, 



lim 7 



J .4Z3 



1 In 2 



(47) 



while 



lim Jrfg (T) 



In- 



d /id 



|G| 



J .4Z3 



2 V 



m .A + "t-b — 1 I In 2 



In 2. 
(48) 



The difference between the mutual information obtained for 
the two order of limits gives 



A-L.R = lim 



Iab(T) ~ T 1I1J i „ 1 AB 



lim 7 (T) 



-im, 



m B — 1 ) In 2. 



(49) 



It is noteworthy that the boundary contribution drops out from 
A7^ B , but that this quantity still has an 0(1) topological con- 
tribution that depends on the total number of disconnected re- 
gions m A + m B of partitions A and B. Hence a topological 
contribution to the entanglement entropy can be filtered out 
directly from a single bipartition using A7^ B , in contrast to 
the constructions in Refs. 5.6 that require a linear combina- 
tion over multiple bipartitions. 



IV. THE TOPOLOGICAL ENTROPY 

Let us now compute the topological entropy using the re- 
sults for the von Neumann entropy in the previous section, 
and the definition given by Levin and Wen 5 



^topo 



lim \—Sh 

r,R — >oc 



1A 
VN 



q2A 



q3A 

"VN 



o4A] 



(50) 



based on the bipartitions shown in Fig. [3] In this specific case, 



= 2, with two distinct £ 



in 



IB 
3B 

From Eq. d35l ), it follows that 



and £ 



and m 2B — 



L 4B 



1, with a single £ 2 p, £3^, £47^ respectively. 
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1 

// 






r , | 






— »■ 



(1) 



(2) 



(3) 



(4) 



FIG. 3: Illustration of the four bipartitions used to compute the topo- 
logical entropy in Ref.|5|. 



2 . . cosh (fy(N-'En, )) 2 smh(fy(N-E-, v )) 

5 topo (r) - sgr = E N 1} M x) 1 L r < ; + E (s l) i^f^ — ^ — ^ 



cosh (^fNj fei 1 cosh {^N 

7 cosh'^^ 



(^iv) V 7 cosh(^iv) 

cosh f4*-(iV- Sap)") , . $mh(^(N-E 3V ) 
(V 3 ) lnx&) ^ ^J- (y^ Iny^ ^ — 



cosh Z£N v ' cosh 



, s cosh(^(N -E 4V )) , x sinhf^(iV-S 47 ,) , 

+ i (4) m* (4) ^ L ^ + (ff« InyW) ^-L , 7 , (51) 

V ' cosh f-^ AM v ' cosh AM 



where 5 t ^ ev is the topological entropy of the GS of the orig- 
inal toric code, which obtains from the — \n(d A d B )/\G\ con- 
tribution to Sy N (T). All the contributions from the (xln x) 
and (y In y) terms cancel since E-^ — E 2 _4 ~ ^3.4 + ^4-4 = ^ 
by construction. 

For finite systems, the known limiting values are recovered: 

S topo (T ^ 0) - S^r - (52) 
S topo (T - oo) - - -In 2. (53) 



In the thermodynamic limit L — > oo and with r, i? kept 
constant, all E a p diverge with the exception of E-^ , which 

corresponds to the inner square of size (R — 2r) 2 in Fig. [3] 
Thus, Eq. (fSTji becomes 



8=1 

- (*( 3 > ln^ 3 > + y<® lny^) e -^ s 3 P + ( £ (4) ln j(4) + g(4) ln ^(4)^ £ -^s 4T , 
= («« ln4 1} + # fctf') + ^E^ - ^E 2 . - ^E 3 „ + ^ 

= (iPlnSP +#ln# )) - ^E 17V (54) 



where we used the fact that 



and that 



N 



J IA ' 
" J 2A ' 



J 3-4 



= e 



4-4, 



E 2 -p 

S 4A> 



E 17 



^4-P' 



(55) 
(56) 
(57) 
(58) 



J 1.4 



J 2A 



J 3.4 



^4A 



(59) 
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to substitute 



- J 2V 



J 4P 



J l-P, 



(60) 



into Eq.d54b. Considering that we are eventually interested 
in taking limit r, R — > oo (i.e., ^> 1), we obtain the 

asymptotic value 



Kitaev L—><x > 
topo 



In 2 



(61) 



for any non-zero value of T, that is thermal equilibrium at any 
infinitesimal temperature leads to a finite loss of topological 
entropy in the thermodynamic limit. In Sec. IVII we discuss 
the implications of this result, and in particular, we propose 
an interpretation that naturally explains why the topological 
entropy reduces to precisely half of its zero-temperature value 
S$f v = 2ba2. 

Notice that Eq. d54l ) is consistent with both the zero- 
temperature and the infinite-temperature limits, Eqs. ( 1521531 ). 
Notice also that the topological entropy in the L — > oo 
limit becomes a pure function of K A Y> 1V ^ j% whose shape 
is illustrated in Fig. [4] The location of the drop, say when 
S topo (T) - S™T = -(^ 2)/2, is given by 



K 



2 ^ 



1 

4' 



(62) 



Even for modest partition sizes with > 100, the drop 

occurs at rather small temperatures and we can approximate 



K A = 
This in turn gives 

e dr °p ~ 



In 



l a.iili ( ^ 
T 



2e" 



T, 



drop 



In 2jE 



(63) 



(64) 



The l.h.s. of the above equation allows for a straightforward 
interpretation in terms of defects in the underlying electric 
loop structure. In fact, e~ x ^l T controls the density of such 
defects in the system, and the equation therefore suggests that 
the drop in topological entropy occurs when the average num- 
ber of defects inside partition \B X becomes of order one. 

In order to understand the behavior of S topo (T) at finite 
temperature and finite system size, notice that the tempera- 
ture parameter K A = — ln[tanh(/3A A )] in Eq. dBTl i always 
appears multiplied by an extensive quantity, be it TV or one 
of the S's. It is therefore convenient to make the reasonable 
assumption that the number of star operators in each subsys- 
tem Ay . . . , A m , and By . . . , B„i B scales linearly with the 
total number of star operators N. Namely, this amounts to in- 
creasing uniformly both L and r, R while keeping their ratios 
fixed, thus simply rescaling the bipartitions in Fig. [3] We can 
then introduce the notation = N^ A , and E p = Nj v , 
with € (0, 1), and j A + £\ j v = 1. Recall- 

ing the definitions of = cosh(K A T, aV /2) and — 




10 10 
K L 12 (arb. units) 




FIG. 4: (Top) Limiting behavior of the entropy difference Eq. l !54t in 
units of In 2 in the thermodynamic limit, as a function of K A T, 1V ^ /2, 
where K A = - ln[tanh(A A /T)] and E 1P ~ (R - 2r) 2 , the area 
of the inner square in Fig. [5] Notice the logarithmic scale on the 
horizontal axis. (Bottom) The same curve represented as a function 
of T/A 4 , for three different values of E 1P = 20, 200, 2000 (from 
right to left). 



smh(K A Yi aV /2), one can replace K A by IC A — K A N and 
all other parameters in Eq. d5"TT l become intensive quantities 
that do not scale with the system size. Temperature and sys- 
tem size are strongly bound together into a single tunable pa- 
rameter K. A in our system. The thermodynamic limit at zero 
temperature is singular, in that the behavior of K. A depends on 
the order of limits. 



A. 



Numerical evaluation of S topo (T) 



The expression for the topological entropy as a function 
of temperature and system size Eq. (IBTt is rather lengthy 
and non-transparent. In this section we illustrate its behav- 
ior graphically, by explicitly evaluating S topo (T) for small 
systems. In Fig. [5] we present the difference (S topo (T) — 



14 



S^f v )/ In 2 as a function of K A 



topo 

tem sizes N 



10 ci ,10 b ,10 9 



K A N, for various sys- 
For convenience, we chose the 




FIG. 5: Topological entropy as a function of K. A for increasing 
system sizes iV = 10 3 , 10 6 , 10 9 . Notice the complete overlap be- 
tween the different curves, due to the fact that the topological entropy 
Eq. d5 1 b becomes a pure function of IC A = K A N when all E's scale 
linearly with N. Notice the logarithmic scale on the horizontal axis. 



values of r and R proportional to V N, so that the above as- 
sumption on the E's holds true, and S topo (T) is a function of 
JC A only. 

In the limit L — > oo, the smooth curves collapse identically 
onto their infinite-temperature value for any non-vanishing 
temperature, and a singularity arises at T = 0. 

The location of the drop is given by K A N ~ 10, from 
which we obtain 



In 



tanh 



A, 



T 



drop 



N 



10 



T, 



drop 



tanh 



(65) 



For large enough system sizes, T drop is small and the above 
equations reduce to 



iVe 



T, 



drop 



In v/TV/5 



(66) 



Once again, the drop occurs when the average number of de- 
fects in the system becomes of order one. (This is consistent 
with the previous result in Eq. (l64l since we made here the 
assumption that all the E's, and therefore E 1P as well, scale 
linearly with N). 



V. THE FULL TEMPERATURE RANGE 

In the regime considered in this paper, finite temperature 
disrupts the cr x -loop structure gradually for finite size systems 



until it is completely destroyed. This happens while the a 1 - 
loop structure is fully preserved, and the topological entropy 
changes overall from 2 In 2 to In 2 (half of the contribution is 
lost). 

The remaining topological entropy should fade away as 
temperature is further increased, and one goes to the regime 
where defects in the erMoop structure also start to appear, for 
a finite energy scale X B . The temperature scale of the drop in 
5 t from In 2 to corresponds to when the distance between 



defects, 



becomes comparable to the system size 



L. (Or equivalently, the average number of defects in the sys- 
tem becomes roughly of order one - compare with Eq. (|64| | 
and (1661) 

It is not obvious how to obtain the exact expression for this 
second step, in contrast with the first step which we calculated 
exactly in this paper within the preserved cMoop limit. Nev- 
ertheless, we believe that the physical picture is the simple one 
(as seen at work in the first drop) that once a handful of defects 
appear in that erMoop structure, the topological entropy will 
plunge much like in the first drop. Pasting the two pictures to- 
gether, we have the two-stage drop of the topological entropy 
sketched in Fig.Q] Clearly, in the limit — X B \ — > the two 
drops are expected to merge together, and in particular in the 
thermodynamic limit the topological entropy entirely vanishes 
for any infinitesimal temperature. 

We would like to point out that a notion of fragility in the 
Kitaev model at finite temperature, in terms of expectation 
values of toric operators, has been discussed by Nussinov and 
Ortiz 16 within their definition of topological quantum order 
(based on gauge-like symmetries). 



VI. CONCLUSIONS 

We calculated the entanglement entropy exactly for the 
toric code at finite temperatures, in a regime where there is 
a broad separation of energy scales between the two cou- 
plings in the problem, X A <C X B . These couplings, from a 
Z2 gauge theory perspective, correspond to the chemical po- 
tentials of electric charges and magnetic monopoles. One can 
define length scales associated with the separation between 



these types of defects, £ 



A,B 



• B ^ T , and for system sizes 
much smaller than the largest of these two length scales, i.e., 
L <C £ B , one of the two loop structures in the system, as- 
sociated with the (T z -basis, is preserved. This is the regime 
where magnetic monopoles are not present in the finite size 
system. In the limit X B — > 00, this holds true for any system 
size. It is in this limit that we obtain the exact result for the 
entanglement entropy as a function ofT/X A . 

Within this hard constrained regime, we find that the en- 
tanglement entropy is a singular function of temperature and 
system size, and that the limit of zero temperature and the 
limit of infinite system size do not commute. The two limits 
differ by a term that does not depend on the size of the bound- 
ary between the partitions of the system into two entangled 
parts, but instead depends on the topology of the bipartition. 
We also calculate the mutual information, obtained from the 
von Neumann entropy by a symmetrization procedure to filter 
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bulk terms at non-zero temperatures and to leave only bound- 
ary and topological contributions. Similarly, the difference 
between the two orders of limits is an 0(1) term that is purely 
topological, depending on the number of disconnected pieces 
of partitions A and B. 

We find that one half of the topological entropy is shaved 
off from its T = value as the temperature increases above 

Tcmls ~ X A / his/N. Above this scale, the loop structure 
associated with the er x -basis is destroyed, while the one asso- 
ciated with the er z -basis survives (recall the X B — > oo). We 
argue that a large but finite value of \ B would introduce an- 

( B) I — 

other scale T cr0 ss ~ A B / In v N, above which the rest of the 
topological entropy should also vanish. 

As these results show, the topological contributions to the 
von Neumann entropy or equivalently to the topological en- 
tropy, are rather fragile for non-zero temperatures. If the ther- 
modynamic limit is taken first, these quantities subside im- 
mediately. However, in practice one should focus on physical 
regimes and not mathematical limits. The reason why these 
quantities are so fragile is that 0(1) defects can destroy them. 
However, one must realize that the length scale associated to 
the defect separation grows exponentially as temperature is 
decreased, and becomes astronomical for temperatures a few 
hundred times smaller than the energy scales A^ B . Hence, 
even if these topological contributions to the entanglement en- 
tropy technically vanish, they are statistically present in large 
but laboratory size physical systems. 

If one is interested in understanding how robust is the topo- 
logical order (information) stored in a single finite system, the 
notion of a statistically non-vanishing topological entropy nat- 
urally translates into the presence of a characteristic time scale 
over which topoogical order is preserved. Such time scale is 
associated with the Boltzman probability for the appearance 
of a defect, namely N e~ XA / T , where Af is the total num- 
ber of degrees of freedom in the system. In sight of a possi- 
ble practical use of such topological quantum information, it 
would therefore be of great importance to compare this persis- 
tence time scale with the one associated to the preparation of 
the system into a topologically ordered state. Preliminary re- 
search in that direction can be found in Ref.[Io|and in Ref.[l7l 

At a more fundamental level, our results suggest a simple 
pictorial interpretation of quantum topological order, at least 
for systems where there is an easy identification of loop struc- 
tures as in the case here studied. Recall that we start from a 
zero-temperature system exhibiting quantum topological or- 
der associated with the presence of two identical underlying 
closed-loop structures. In particular, the corresponding topo- 
logical entropy equals InD 2 , where D = 2 is the so-called 
quantum dimension of the system. By allowing one of the 
two loop structures to be thermally disrupted, and by raising 



T — > oo while the other loop structure is fully preserved, we 
arrive at a classical system with a single (therefore classical) 
underlying loop structure, and exhibiting precisely half of the 
original topological entropy (In D). This is strongly sugges- 
tive that (i) the two loops structures contribute equally and 
independently to the topological order at zero temperature; 
(ii) each loop structure per se is a classical (non-local) ob- 
ject carrying a contribution of In D to the topological entropy; 
and (iii) the quantum nature of the zero-temperature system 
resides in the fact that two independent loop structures are 
allowed to be superimposed and thus coexist in the system. 
In this sense, our results lead to an interpretation of quan- 
tum topological order, at least for systems with simple loop 
or membrane structures, as the quantum mechanical version 
of a classical topological order (given by each individual loop 
structure). 

Finally, we would like to comment on the fact that the same 
0(1) defects that deteriorate the topological entropy of the 
system should also deteriorate its usefulness for topological 
quantum computing. A handful of stray unaccounted defects 
winding and braiding around others that are accounted for 
in the computational scheme will lead to errors. These de- 
fects can be thermally suppressed, if the temperature is small 
enough and the system not too large, so that unwanted defects 
have a small probability of appearing in the sample. Thus, 
quantifying topological entropy at finite temperature and fi- 
nite system size is meaningful in quantifying, in a statistical 
sense, the degree with which a physical (finite) system retains 
topological order. 

Although the results presented here were derived in the case 
of one of the coupling constants being infinite, we have re- 
cently been able to extend the calculations to the case where 
both coupling constants are finite [C. Castelnovo and C. Cha- 
mon, in preparation]. The two contributions to the topologi- 
cal entropy due to the underlying gauge structures are shown 
to behave additively, and indeed the behavior conjectured in 
Fig- His confirmed. 
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